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ABSTRACT 


A stochastic stellar atmosphere model which simulates the behavior of pho- 
tons in the atmospheres of stars by following individual photons is described and 
tested. The model accounts for stellar turbulence, atmospheric curvature, stellar 
rotation, collisional excitation, collisional de-excitation, and stellar winds. Testing 
shows that this new method of stellar radiative transfer is successful in simulating 
all of these processes. Particular results of interest include simulation of solar limb 
darkening and spectral synthesis of three different solar absorption lines; Ti 
6258.1 A, K 4044.1 A, and Sr 14607.3 A. 
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I, INTRODUCTION 


Recent advances in spectroscopy have led to a heightened interest in radiative 
transfer theory. As R P. Kudritzk stated (Ref. 1) 


sive space experiments or extensive ground anes observations using 


oa he foci and enggreingr on ow ings aw eeveopaen 

What is lacking is a comprehensive, workable, quantitative, radiative transfer 
theory with which to extract as much information as possible from available stellar 
spectra. The principle problem faced by modern radiative transfer theory is the 
treatment of non-local thermodynamic equilibrium (NLTE) portions of stellar atmo- 
spheres. A few examples are: B- and A-type supergiant stars where NLTE theory 
is needed in the photosphere and the wind region, solar-type stars where strong 
absorption lines occur in the upper photosphere and the chromosphere; pre-main 
sequence stars with high mass loss and disks that protrude into the photosphere; 
and, probably the most complex case, O-type stars, where all LTE assumptions 
breakdown completely even deep in the photosphere. 

Many different theoretical techniques have been applied to the problem of 
radiative transfer with varying degrees of success. In this paper, a new approach 
which follows individual photons through the stellar atmosphere of the star is 
examined. 





A. BACKGROUND ON RADIATIVE TRANSFER THEORY 


Radiative transfer theory is the study of how photons propagate through a 
medium. The theory has its origins in astronomy and the study of star atmo- 
spheres and spectra. But there are many additional military and civilian applica- 
tions of radiative transfer theory (i.e., fusion bomb triggering, nuclear radiation 


transfer, atmospheric sciences). 


To describe the interaction of matter and light it is necessary to characterize 
the state of the matter as well as that of the radiation field. Typically, the 
assumption of thermodynamic equilibrium (TE) or loca] thermodynamic equilibrium 
(LTE) is made to simplify this characterization. These assumptions significantly 
reduce the complexity of many of the radiative transfer computations. A quick 
review of these assumptions is in order. TE is the most restrictive case. A gas in 


TE will satisfy the following conditions: 


1. Each species of gas particle will have a Maxwellian velocity 
distribution specified only by its mass and a temperature; 

2. The distribution of the excited states of each ion will be 
accurately described by Boltzman’s law specified only by the 
energies of excitation and a temperature; 

3. The distribution of the ionized states of the atoms is accur- 
ately described by Saha’s equation, again specified only by the 
energies of ionization and a temperature; and 

4. Radiation intensity is accurately described by Planck’s function 

specified only by a temperature. 


It is important to note that in TE, the same temperature variable is utilized 


in all four relationships given above. 





TE can only be achieved in an isolated system, so that the gas can not inter- 
change any energy or mass with its surroundings. It is obvious that this condition 
is rather stringent and difficult to meet in a stellar atmosphere. 

Since true TE cannot be fully justified, the less stringent LTE is usually 
invoked. In this second case, the gas wil] have TE characteristics over a loca! tem- 
perature zone. This assumption is justified when collisions among atoms of the gas 
occur rapidly enough to be the dominate cause of excitation. When this occurs, the 
kinetic temperature of the atoms in the local region can be used to accurately 
describe the gases’ atomic properties. The fourth assumption will not hold for the 
LTE case because radiation escaping from deeper regions of the star cannot be 
described by the kinetic temperature of the local region; therefore, the radiation 
field can deviate arbitrarily from the Planck distribution. 

The LTE assumption is a good description of the state of the matter deep in 
the photosphere of a star, where the majority of the energy is transfern-d via con- 
duction or convection. This is for two reasons: 1) due to the high opacity, little 
radiation (energy) can escape directly, thus preserving energy conservation; and, 2) 
use of Boltzman’s equation is justified because collisions are the main mechanism 
for atomic excitation thus allowing the excitation to be described by a single tem- 
perature. The LTE assumption begins to break down as the collisional transition 
rates of the ions approach the radiative transition .ates. This decrease in colli- 
sional rates relative to the radiative rates occurs as the photons work their way up 
through the atmosphere during which the density of the gas and the kinetic tem- 
perature of the atoms decreases as energy is radiated away. 

When the LTE assumption is no longer valid, the state of the gas can be cal- 
culated by statistical equilibrium (SE) methods. This approach balances the 
number of transitions due to all processes into and out of a given atomic level. This 
method is very difficult given the large number of energy levels and ionization 
states available to an atom. The method of statistical equilibrium is always valid 














and the assumptions of LTE and TE are only special cases of SE. 


With an understanding of the characterization of the medium, the radiative 
transfer equation is 


dl, 
dw 





=I,-S, (1) 


where: /, = initia) intensity as a function of frequency, 
S, = source function as a function of frequency, and 
t, = optical depth as a function of frequency. 

t, and S, are given by 


L 
= [mpde , and (2) 
Jy 
Sy,=— 3 
a (3) 


where: dx = geometrical depth, 
L = path length, 
p = density, 
2 
x, = extinction coefficient ( £~ ), and 


Jy 2 emission coefficient ( ES ). 
g rad 


The simplicity of equation (1) is deceiving. The complexity of the radiative 
transfer problem lies in determining the properties of x, and j, . «, and j, are 
dependent on the intensity which renders equation (1) a first order non-linear 
differential equation vastly complicating the problem. 

Radiative transfer theory has been successful in explaining many features of 
stellar spectra. There are some important pmblems in which the restrictive 
assumptions of this theory fail to adequately describe the stellar atmosphere. 








1. NLTE 


Some of the more important spectral lines (i.e., strong absorption lines 
due to H, He, and metals) have non-LTE (NLTE) characteristics. This is because 
the line centers are formed high in the atmosphere where the density and kinetic 
temperature are low and radiative processes dominate over collisional processes. 
This NLTE behavior is not restricted to exotic high-temperature, high-luminosity 
stars; the sun contains several NLTE lines of interest. The importance of under- 
standing NLTE line formation can not be over emphasized. Spectral measure- 
ments of faint stars +-- ‘cally can onlv detect strong absorption lines over back- 
ground noise. Many of these absorption lines are NLTE in nature. Therefore an 


understanding of this NLTE line formation region is of paramount concern. 


2. Winds 


Radiatively driven winds in stellar atmospheres are an additional prob- 
lem for the LTE assumptions. Many stars, particularly high-temperature, high- 
luminosity, very young, and very old stars, have rapid mass loss (typicelly 1 earth 
mass/yc) via radiatively-driven winds. These winds produce expanding stellar 
atmospheres where the LTE assumptions quickly break down. These expanding 
atmospheres also preferentially Doppler shift line profiles, which complicates the 
radiative transfer problem. Loss of matter must also be considered when analyzing 
the validity of the LTE assumption. Boltzman’s equation requires a closed system 
in which the total mass of the system is conserved. Rapid mass loss is in violation 


of this requirement. 


B. CURRENT STELLAR ATMOSPHERE MODELING 


Stellar atmosphere modeling is an attempt to reconstruct the structure of the 
stellar atmosphere based on spectral observations and physical laws. The complex- 


ities of the dynamics of the atmosphere as well as the complexities of radiative 





transfer theory prohibit direct measurement of many of the physical parameters. 
Therefore, a model is hypothesized to generate the observed spectral information. 
Once the model has been constructed which reproduces the observed information, 
additional observations can be made to further test the model. While this may not 
produce unique atmospheric parameters, stellar atmosphere modeling has become 
an indispensable tool in the analysis of stellar spectra. Most models contain the 


following simplifications: 


1. LTE is assumed throughout the atmosphere; 

2. Plane-parallel geometry is assumed ignoring atmospheric curvature; 
3. Hydrostatic equilibrium is assumed neglecting any mass losses; 

4. Fine structure is negligible (i.e., granulations, flares, sunspots); and 
5. Magnetic fields are excluded. 


Stellar modeling can be divided into three distinct types: continuum, line, 
and spectral synthesis. Each method will be examined further in the next three 


sections. 


1. Continuum modeling 


Continuum modeling seeks to reproduce the characteristics of a stars 
continuum radiation. The source function is calculated from the Planck distribu- 
tion at a given depth of formation. This quantity is referred to as the surface flux. 
The extinction term from the radiative transfer equation (eq 1) is based on contin- 
uum absorption coefficient calculations and laboratory measurements of atomic 
properties (ie., H , H™ , electron scattering). The standard radiative transfer 
integral is then solved, producing the calculated spectrum. The success of these 
types of models is remarkable considering their simplistic assumptions. Figure 1 


shows a fit according to Melboure [Ref. 2). 


Continuum model comparison to Observations 


51 Peg 
o Continuum 
» Observed 
— Mode! 





Figure 1: The observed energy distribution in 51 Peg(G5V) corrected for line 
absorption compared with model calculations [Ref. 2]. 


Note that the observed data has been adjusted to account for line absorption at the 


shorter wavelengths. 


The continuum modeling technique contains three major problem areas: 
1) absolute continuum flux calibration (i.e., where is the continuum located), 2) 
NLTE effects on the temperature gradient, and 3) interstellar reddening. The first 
two of these problems is intrinsic to the continuum modeling technique and not 
easily solved; interstellar reddening complicates comparisons to observed contin- 


uum. 


2. Line modeling 


Line modeling efforts have been directed at determining the shape of line 
absorption profiles. With accurate profiles much can be learned about stellar atmo- 
spheres. The line absorption profile is dependent on many important atmospheric 


parameters (i.e., pressure, temperature, radiation field, rotation, winds, turbulence, 
depth of line formation, etc.). It is the interaction of these parameters that pro- 
duces the absorption profile and thus reveals information about the conditions at 
formation of the line. In addition, line modeling provides information about the 
processes that produce various line profiles. 

The radiative transfer equation (eq 1) is the basis of these calculations. 
But the central problem remains how to obtain an accurate source function. Typi- 
cally, S, is approximated by a combination of a line-source function and a 
continuum-source function. The line-source function is calculated from LTE con- 
siderations and the continuum-source function is calculated from surface flux calcu- 
lations. These calculations have produced good results when compared to actual 
absorption-line profiles. Figure 2 and 3 show the results of line modeling done by 
Carbon and Gingerich [Ref. 3), and Curtis and Jefferies (Ref. 4) respectively. 


Line models compared to observations 


0 Gray and Evans (1973) 
» Wright et al. (1964) 
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Figure 2: The observed H,, profile in 9 Leo compared with LTE line profile cal- 
culations. The model calculations fail to correctly model the line cores (Ref. 
3). 








Figure 3: The observed Na D line profile compared with LTE line profile cal- 
culations. Again the mode) calculations fail to correctly model the line cores 
(Ref. 4). 


The depths of the observed profiles are deeper than those calculated by 
line-modeling theory. This is one of the inadequacies of the current line-modeling 
theory. Gray cites the effects of turbulence and NLTE population levels as typi- 
cally deepening line cores (Ref. 5). Inclusion of these phenomena improves line fits 
markedly but some differences still remain. The line-modeling-technique limita- 
tions are not only confined to NLTE effects at line cores. Weak-line calculations 
also fail to match theoretical curves of growth. Curves of growth describe the 
theoretical growth of an absorption line as a function of abundance, temperature, 
and pressure. Weak-line calculations fail to match these curves although inclusion 
of turbulence improves the fits. 


3. Spectral synthesis 
Spectral synthesis is a method of spectral analysis developed to calcu- 
late the chemical abundance in stellar atmospheres. High-resolution observations 
of spectral lines are compared to the synthesized lines. The individual abundances 
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of the atoms are adjusted until a good fit in the synthesized line is achieved, thus 
determining the chemical composition of the star. The intensities of the spectral 
lines are calculated by evaluating the standard radiative equation (eq 1). The 
results of some spectral synthesis are shown in Figure 4 and 5. 


Spectral synthesis compared with observations 


Computed 


0.8 


06 
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Figure 4: An observed portion of the solar spectrum compared with Ross and 
Aller spectra] synthesis [Ref. 6). 


Figure 4 is a comparison of computer synthesized lines to observed solar lines [Ref. 
6]. The fit is fairly good, although turbulence effects have been added. Figure 5 is 
the synthesized spectrum of a carbon star (TX Psc) compared to observations [Ref. 
7). Although the fit shows significant deviation from observations, Figure 5 shows 
how even very crowded regions of the spectrum can be synthesized. 

This third approach contains many of the same problems as the line- 


modeling method. The emphasis in this method has been on the problem of line 
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Figure 5: The observed spectrum of TX Psc compared with Synthesized spec- 
trum [Ref. 7). 


blending in crowded regions and not on the line profile itself. 


4. Limitations 


Although these modeling methods have been very successful, each does 
possess some limitations as noted. In addition, all of these methods contain four 
generic problems: 1) the assumption of a plane-parallel geometry, inappropriate for 
large extended stellar atmospheres; 2) they lack an easy method by which to incor- 
porate stellar winds and rotation; 3) Line profiles and continuum models fail to 


reproduce this observed limb darkening behavior* across the solar disk; and 4) 
LTE atmospheres are assumed. 


* Observations of the sun’s disk reveal that intensity decreases as observations near the 
solar limb. 











5. Other approaches 


In an effort to solve the NLTE problem, Kudritzki (Ref. 8] used the 
method of linearization proposed by Auer and Mihalas [Ref. 9] to develop a com- 
puter code designed to model the structure of the NLTE atmosphere. The method 
accounts for all possible coupling between variables at all different depths within 
the constraints of statistical, hydrostatic, and radiative equilibrium. This is essen- 
tially the statistical equilibrium method. To show the complexities of these calcula- 
tions, the O II Grotrian diagram used by Becker and Butler [Ref. 10, 11] is shown 
in Figure 6. 


Grotrian diagram 





Figure 6: O II Grotrian diagram showing the 200 different transition which 
must be calculated [Ref. 10, 11). 


The computer code must be given, or calculate, all atomic parameters and calculate 
all transition rates into and out of these 200 different atomic energy levels. 
Kudritzki (Ref. 12] has demonstrated that a large number of transitions are neces- 
sary for acceptable accuracy. With accurate atomic parameters, the results of this 
method can by highly accurate (Figure 7). Discrepancies in the line core were 
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corrected for by Herrero [Ref. 13] by including new collisional broadening terms 
and coupling between H and He II. 


NLTE line profile compared to observations 
1 1. 
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Figure 7: Observed H, profile of the O 9.5 V star t Sco compared to NLTE 
calculations {Ref. 9] (left) and recent improvements by Herrero [Ref. 13] 
(right). 


In an effort to explain multi-line wind transfer problems, Abbott and 
Lucy [Ref. 14] developed a Monte Carlo technique for the exploration of stellar 
wind effects on the spectra from stars with strong stellar winds. The atmospheric 
parameters used by Abbott and Lucy were from Kudritzki, Simon, and Hamann 
(Ref. 15]. Their technique follows energy packets, which may be as small as hv, 
from the top of the photosphere to their emergence through the stellar wind por- 
tion of the upper atmosphere. In this way wind effects can be simulated and 
wind-driving forces explored. Results of this technique are shown in Figure 8 and 
reveals general agreement with observations of € Pup in the ultraviolet. 
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Wind simulation compared to observations 
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Figure 8: Abbott and Lucy fit to observed ¢ Pup spectrum. Prominent features 
are identified on the top axis [Ref. 14]. 


C. STOCHASTIC PHOTON-FOLLOWING SIMULATION 


1. Motivation for SAM 
In view of the limitations with current atmospheric modeling and radia- 
tive transfer methods a new approach has been explored by W. B. Weaver at the 
Monterey Institute for Research in Astronomy. He has developed a stochastic 
Monte Carlo stellar atmosphere model (SAM). The advantages of this approach 
are listed below: 


1. Interaction of photons and matter can be modeled explicitly and in detail 
permitting easy inclusion of NLTE effects; 
2. Computations are done in spherical geometry; 


3. Influence of stellar rotation is included; 
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4. Influence of stellar winds is included; 


5. 2 4 dimensional! geometry is used naturally reproducing limb darkening; and 


6. Output is in 2 dimensions; 4 and 6 ; where 0 is the angle from the surface 


normal. 


SAM is a true simulation, with the fundamental atomic processes being 
modeled. The advantage of this approach allows the modeler to investigate funda- 
mental physical processes, leaving the radiative transfer computations to the com- 
puter. 


The theoretical concept behind SAM is simple. At large optical depth, the 
radiation field is accurately described by the Planck function. Individual photons 
selected from the Planck distribution are followed as they migrate through the 
atmosphere. The properties of the photosphere are then accurately described by 
the LTE assumptions except, as stated before, for the radiation field. But by keep- 
ing track of the photons from the deeper region of the atmosphere, where the 
Planck function is valid, an accurate picture of the radiation field can be produced. 
In other words, SAM is able to reproduce both the radiation field as well as the 
state of the atoms in the important LTE region. 

If one then continues into the NLTE region, the radiation field deter- 
mined in LTE will help to formulate the characteristics of the NLTE portion of the 
atmosphere, thus approximating a solution to the NLTE radiation transfer prob- 
lem. 

Stellar wind effects can also be simulated. The outer atmosphere of SAM 
can easily be set in motion. This allows modelers to gain insight into the effects of 


rapidly moving atmospheres on line profiles. 
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2. SAM1_13 


SAM1_13, the most recent version of SAM, is used as the starting point 
for this thesis. The method employed by SAM1_13 to simulate the radiative 
transfer problem is to generate a number (usually 500,000 photons/A) of photons 
and distribute these photons evenly over the wavelengths in question. The 
assumption being the Planck function is flat over a few angstroms. 


Photons are started at a calculated thermalization depth, usually deep 
within the photosphere. The calculations of thermalization depth are based on con- 
tinuum photons with optical depth of 10. The individual photons are now permitted 
to migrate through the atmosphere. Random absorption probabilities are produced 
for continuum and line processes. Since continuous variation of atmospheric pro- 
perties is impractical, atmospheric parameters are held constant for one to ten per- 
cent absorption probability. These small absorption probabilities are added until 
the initial randomly chosen absorption probability is reached. This procedure is 
followed for both line and continuum processes. The lengths traveled for the two 
processes are then compared, with the shortest being chosen by the computer as 
the absorption process the photon undergoes. At this point the photon is con- 
sidered absorbed and a re-emitted photon is produced. For simplicity, collisional 
processes are not considered. The new photon is redistributed in wavelength and 
emission angle according to the absorption process involved (i.e., H~ absorption, 
line absorption). This process is continued for the photon until an assumed escape 
altitude is reached, usually the top of the photosphere, or it is destroyed by pro- 
pagating back into the thermalization depth. Radiated photons are then collected 
into bins chosen by the user, typically of width .01 angstroms. Angle to the surface 
normal at the time of departure is also retained and radiated photons are binned 
accordingly. 

SAM1_13 requires two programs to be run prior to actual photon genera- 


tion. These programs read in atmospheric parameters and generate constants that 


16 











are needed by SAM. Figure 9 shows the order in which the different programs are 
run. Figure 10 shows the conceptual flow diagram of SAM’s logic. Figure 11 is a 
detailed diagram of the photon absorption loops as described above. 


SAM PROGRAM SEQUENCE 













OWEN: Takes atmospheric parame- 
ters (i.e., continuous opacity, tempera- 
ture, density, ...) and arranges them 
into the correct format 


NSUM: Generates spline fits for all 
parameters received from OWEN. 
Also calculates Voigt a, and collisional 
rates and stores all these parameter 


SAM: Propagates photons through the 
generated atmosphere. 
VEL: Applies rotational Doppler 
broadening to SAM output. 






Figure 9: Order of execution of programs. 
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SAM CONCEPTUAL FLOW DIAGRAM 


Input spectral region of interest and 
atomic parameters 


Begin calculation cycle (blackbody pho- 
tons + collisional photons) 


Start blackbody photons 






a) Calculate propagation distance for 
line and continuum processes 


b) Does photon escape? 
c) Does photon re-thermalize? has 


aaa 
Is blackbody cycle complete? 


yes 







a 


Figure 10: The logic flow of SAM during radiative transfer. 








LINE ABSORPTION LOOP CONTINUUM ABSORPTION LOOP 
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Figure 11: The logic flow of the two absorption loops in SAM. 


Spherical geomet.y is the natural choice for photons propagating in SAM. 
The position of photon’s can be described by the three spherical cc ordinates p, 8, , 


and 6,. Stellar atmospheres have spherical symmetry rendering p the only impor- 


tant coordinate in describing the photon’s position. The photon’s direction of pro- 


pagation can be described by the two coordinates 6 and @. Where 6 is the angle 


from the surface normal (p), and @ is the remaining angle defining the orientation 


of tne plane formed by 6 and p. The geometry is simplified by assuming that the 


orientation of this plane ($ angle) is symmetric for all photons. The @ dependence 
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can be condensed into the random @ calculations. With the condensation of 6, the 
coordinate system becomes a floating reference frame. 
Geometry of SAM 





Figure 12: Orientation of p, 8, and ¢@. 


In order to minimize computation time, all photons being radiated are 
assumed to propagate in the direction of the observer. At first, this may seem to 
be an inappropriate assumption, however, once the photons have left the atmo- 
sphere of the star, the © dependence is removed and photons may be redistributed 
to propagate in the direction of the observer. Another way to look at this 
simplification is to think of the photons traveling backwards from the observer. 
When the photons are first absorbed at the upper surface of the star the incoming 
® value is completely lost. It is at this point where ® can have any value, therefore 


a 6 is chosen so the photon propagates in the direction of the observer. 
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8. The shortcomings of SAM1_13 


Although SAM1_13 generates some interesting results, it has several 
shortcomings. 1) It fails to account for an unequal number of photons being pro- 
duced and destroyed in the atmosphere above the thermalization depth by colli- 
sional processes. 2) It fails to account for energy conservation in the upper photo- 
sphere. This failure caused the upper atmosphere atoms to have LTE population 
levels when NLTE population levels would be more accurate. 3) It fails to have a 
mechanism for imposing stellar rotation effects onto the radiated photons. Rota- 
tion effects can be significant in some stars (i.e., line broadening of up to several 
angstroms). And 4) it contains simplified wavelength redistribution functions. In 
addition, it is felt that due to the large quantity of photons being generated and 
the many interactions each photon undergoes, a more rapid method of determining 


current atmospheric parameters is needed. 
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II. OBJECTIVES OF THIS THESIS 


The main objective of this thesis is to improve SAM by including the following 
physical processes to provide a more complete model, including some NLTE effects. 
The physical processes being added are: stellar rotation, collisional excitation, colli- 
sional de-excitation, optimized photon starting depth, and an improved photon 
redistribution function. Additionally, minimizing the time required to calculate stel- 
lar atmosphere parameters will also be added. Once the improvements are made, 
SAM-generated lines are tested against three observed solar lines. These lines are 
Ti 16258.1 A, 3d°4s — 3d54p, multiplet # 104, J 4-5; K 14044.1 A, 4s-5p, multiplet 


#3, d + - 13; Sr 14607.3 A, 5s* — 5s 5p, multiplet # 2, J 0-1. In addition, a com- 


parison of the SAM-generated limb-darkening with the observed solar limb- 
darkening is made. Finally, an attempt is made to reproduce the characteristic 
structure observed in the spectra of stars with strong non-thermal winds. 
Although NLTE features are programmed into SAM (i.e., two temperatures for the 
outer photosphere), no attempt is made in this work to generate NLTE outer atmo- 
sphere atomic population levels. 











Ill. IMPROVEMENTS TO SIMULATION 


A. STELLAR ROTATION 


Doppler shifting due to stellar rotation can significantly broaden observed line 
profiles. Since this is due to the integration of shifted, unbroadened lines from all 
portions of the disk visible to the observer, the total amount of light removed from 
the continuum is unchanged. At high rotational velocities, a few hundred kilome- 
ters per second, additional line broadening can be as much as several angstroms. 
Inclusion in SAM is necessary to correctly reproduce line profiles from stars rotat- 
ing faster then 1 km/sec. The star is assumed to rotate as a solid body. This 
assumption is justified for many stars. For example the differential rotation for 
the sun from the equator to the 75° latitude is only 25%. Because SAM follows 
individual photons as they travel through the rotating stellar atmosphere, Doppler 
shifts within the atmosphere (internal Doppler shifts) must be examined as well as 
those affecting the observations of an external observer (external Doppler shifts). 

There are no internal Doppler shifts because of the solid body assumption. 
Since two internal points in the atmosphere have no relative motion, no Doppler 
shifting can occur to a photon traveling between these two internal points. 

The only broadening due to rotation comes from external Doppler shifts. Fig- 
ure 13 shows photons escaping from the limbs of the star which will appear to the 
observer as blue and red shifted with respect to the center of the disk. 


Each exiting photon in SAM must be Doppler shifted according to its exit 
location. The output of SAM is divided into angular increments. A rotating body 
inclined at an angle i to the line of sight of the observer can have its surface 
divided into a series of strips of unique radial velocities. Superimposed over the 
angular divisions of SAM are these unique radial velocity strips. In a separate 
software routine called VEL, the geometrical regions are calculated and the indivi- 
dual photons from each region are then Doppler shifted as appropriate. The 
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Photon Doppler shift diagram 
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Figure 13: Photons are Doppler shifted as seen by the observer depending on 
the portion of the star from which they radiate. 


general equation for this procedure is 


we Tr UF x 
= ——— 4 
NP\v,r ) ee A aes Pre @)?’ (4) 


where NP(v,r) = number of photons as a function of v and r, 
A, = ratio of total area to partial area of rings and velocity strips, and 
Uri = ratio of strip velocity to speed of light. 


The input rotational speed for VEL is equal to v sin(i) . Figure 14 shows the 
geometric regions of interest and the complicated areas that must be calculated. 





Velocity and angular regions 


onguicr 
increments 


velocity 





apparent disk of the stor 
Figure 14: Velocity and angular regions as projected onto the surface of the 
star. The regions that must be calculated are defined by the intersection of 


the these geometric regions. 


Figure 15 shows the results of applying VEL to Ti 46258.1A for three different 
rotational velocities. These results agree with theory and show that rotation 
effects can be quite significant. Solar rotation is 2 km/sec, and has little effect on 
line profiles. 
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Stellar rotational broadening 
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Figure 15: Different rotational velocities as applied to Ti 46258.1A The top 
line is for a rotational velocity of 50 km/sec, the middle line is for rotational 
velocity of 10 km/sec, and the bottom line is for 1.0 km/sec. The flux offset is 


induced for viewing clarity. 


B. COLLISIONAL EXCITATION AND DE-EXCITATION 

A major step towards making SAM a fully NLTE simulation is the introduc- 
tion of collisional processes. Photons need to be added as a result of collisional 
excitation of an atom followed by photo-de-excitation. In addition, photons need to 
be removed as a result of photo-excited atoms which undergo collisional de- 
excitation. 

The key to each of these processes is the value of the collisional excitation 
(C;,) and de-excitation (C,,) rates. From Jefferies (Ref. 16], a semi-empirical 
equation based on the dipole approximation provides a means to calculate the exci- 
tation rate (C;,,) from the lower to the upper state. 
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Cy, = 2.1607! exp (-0) Tn, (5) 


— Jo 
where a= KT ’ 


No = energy of excitation, 

n, = electron density, 

f = oscillator strength, and 
T = temperature. 

With the value of C;, known, calculation of C,, follows from the principle of 
detailed balance. This principle states that, in TE, all processes occur exactly as 
rapidly as their inverse. Although TE is required to derive this equation (eq 6), the 
fact that it only contains various atomic parameters that are properties only of the 
atoms themselves allows the equation to be valid under NLTE conditions as well. 
The equation 


= Foye 
Cu : exp( KT Cru (6) 
where g; = statistical weight lower level, and 


& = statistical weight upper level 


can be used to relate these two collisional rates. 


To compute the number and starting depth of the photons produced by the 
collisional excitation process, the sum of all the photons produced by this collisional 
process is first calculated from the starting depth of the blackbody photons (z9) to 
the exit altitude of the atmosphere (z,). The equation for the total number of colli- 
sional photons from depth 29 to 2; is 


Ne, = ECyy (dz (7) 
zo 
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and for the blackbody photons is 


Nip = 2ne. 1 (8) 





where: Nc , = Total number of collisional photon /cm? from zo to zr , and 
N,, = Number of blackbody photon / cm? . 


In this way a ratio of the number of blackbody photons to collisionally- 
produced photons can be found. For example, if the ratio of blackbody photons to 
collisionally-produced photons is 1 to 3, for every blackbody photon SAM generates, 
it must also produce three collisionally-produced photons. A collisional photon vs. 
depth distribution is also calculated during the summing of the collisional photons. 
Collisionally-produced photons are assigned a starting depth randomly based on 
this probability distribution. 

Photon removal caused by collisional de-excitation is calculated in the line 
interaction subroutine. If a photon is absorbed by an atom, then a collisional de- 
excitation probability is checked against a random number and those photons that 
fall within the de-excitation probability are removed. 

In this way, the collisional source function, dependent on depth, and its 
accompanying extinction process are simulated. Unlike typical treatments by radi- 
ative transfer theory, SAM’s source function does not need to include photo- 
excitation and de-excitation. These processes are inherent by the way SAM follows 
individual photons. 

Figures 16 and 17 show the SAM-generated Ti (6258.1 A line with and 
without the collisional source function included. Note the line without the colli- 
sional source function (Figure 16) is deeper than the line with the source function 
included. This effect is caused by an imbalance between photons being produced 
by collisional excitation and photons being destroyed by collisional de-excitation. In 
the outer photosphere there are more photons added by collisional excitations than 
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removed by collisional de-excitation. 
Ti line with no collisionally excited photon 


flux 


m—u No collisional photons 
6259.0 
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Figure 16: SAM generated line with no collisionally-produced photons added. 
Ti line with collisionally excited photons 
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6259.0 


m—u Collisional photons odded 
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Figure 17: SAM generated line with collisionally-produced photons added. 
Note that the line depth has decreased due to the increased number of photon 





from collisional excitation. 
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C. INTERPOLATION TABLES 


In an effort to reduce the program’s computation time, an analysis of SAM’s 
computation properties was undertaken. The results indicated the computer, a 
Sun SPRAC1 workstation, was spending a significant amount of CPU time calcu- 
lating and re-calculating atmospheric parameters. Specifically, the subroutine 
STAR, where spline fits to atmospheric parameters were being calculated, was the 
largest user of CPU time. 

A simple re-design of the subroutine, replacing the spline fits with interpola- 
tion tables, was completed. The interpolation tables are calculated in NSUM using 
the same spline fit routine. The interpolation tables are incremented in 1 km dis- 
tances. This distance may be varied when stars other then the sun are examined. 
To ensure no detail in the structure of the atmosphere is lost, graphs of all atmos- 
pheric parameter were compared to the interpolation tables (see Appendix A). The 
net effect of this re-design is a 300% increased in algorithm efficiency with negligi- 


ble decrease in accuracy. 


D. OPTIMIZATION OF PHOTON STARTING DEPTH 


Another significant degradation to computation speed is the chosen photon 
starting depth. It was found that photons which started too deep within the star 
would undergoes needless absorption and re-emission. A starting depth which 
minimized computer time is in order. The depth chosen is where the line absorp- 
tion and continuum absorption processes had identical optical path lengths. Above 
this point in the atmosphere, the difference between these two optical path lengths 
shaped the outgoing photon flux. Below this point, the differences in the path 
length may shape the outgoing flux, but the flux is completely re-ordered as it 
migrates through the equal path length region. Figure 18 shows Ti (A6258.1 A ) 
line generated at the optimal depth (484km) and one generated with z equal to 
490km. 
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There is little different between the two line profiles, but the computer run time 
for the line profile with z equal to 484km is 10% as fast. 


Ti line generated at two different starting depths 
3500.0 


flux 





@-© storting depth of 22490 km 
Oo starting depth ot 7=4184 km 


0.0 pe teee ee ee 
6257.0 6258.0 6259.0 


wovelength 
Figure 18; The same Ti line started at two different starting depth. Note that 
the line shape is almost identical. Difference in flux of the two lines is in- 


duced for viewing clearity. 


E. REDISTRIBUTION 


Photons undergo many different absorption and scattering processes while 
propagating through the stellar atmosphere. How these photons are redistributed 
in frequency is of great importance in determining the final line profile. SAM1_13 
had two simplified redistribution functions; complete Compton-scattering redistri- 
bution for the continuum absorption process, and modified complete Doppler redis- 
tribution for the line absorption process. These two simplification are unsuitable 
for accurate line profile simulation. Improvement to both frequency redistribution 


functions is necessary. 
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The major source of continuum absorption in the sun is from H~ . H™ has a 
low ionization energy at .754 ev, therefore photons with wavelengths less then 
16,450 A have enough energy to ionize H~ . Compared to other sources in the visi- 
ble region of the continuum, H™ is stronger by a factor of 4 than that due to neu- 
tral H and much larger than that for other processes. Therefore, for the sun, all 
continuum absorption is assumed to be caused by H™ . This assumption is valid 
for many types of stars. H™ low ionization energy results in a simple redistribution 
function. All photons absorbed by the continuum process are completely redistri- 
buted in frequency. In consideration of minimizing computation efforts, photons 
redistributed out of the wavelength region of interest are replaced by the introduc- 
tion of a random wavelength photon from these other regions. In this way, photon 
flux is conserved. 

The line-absorption frequency-redistribution function is a combination of 
Doppler redistribution (Gaussian in shape) and collisional redistribution 
(Lorentzian in shape). It is the combination of these two distributions that result 
in the familiar Voigt profile for line emission and absorption. The rate of collisions, 
the Voigt a parameter, the kinetic temperature, and the spontaneous emission 
coefficient (Einstein A), are all important factors in determining line-absorption 
frequency-redistribution. Doppler redistribution is based on atoms having Maxwel- 
lian velocity distribution [Ref. 17]. The standard Gaussian profile is given by 


~«42g7 
P(x)=* 


ao’ (9) 





where: o = half width at half maximum. 


The Doppler redistribution function from Hummer [Ref. 18] is given by 


iii et 4d 


. 1 
Gi (10) 





If we let Ay = V2o, then 
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where io = Ay = “47 


Xo = initial wavelength, 
c = speed of light, 
K = Boltzman constant, 


= P(x) (11) 





T = kinetic temperature, 
M = mass of the atom, 
Ag = average Doppler shift. 
For the collisional redistribution, the standard Lorentzian profile is given by 





P(x) = = (12) 
xo(1+—— 
o? ) 
The collisional redistribution function from Hummer [Ref. 18] is given by 
(A"+A° ) 
v) = —_>————_; » (13) 
a R{(Av)?4+{A" +A° )*) 
where A” = average natural broadening, 
A° = average collisional broadening. 
Let o = A" +A° and rearranging 
ov) = + =P(x). (14) 
no(l+=5) 
The Voigt a parameter is 
A" +A° 
a = , 
Ag 
therefore o=Aga. 
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Normal random number distribution was modified by the Box-Muller method 
(Ref. 19] to produce a temperature-dependent Gaussian (eq 15). Random 
Lorentzian distribution dependent on the Voigt a parameter and temperature was 
produced by inverse transform method (eq 16)[Ref. 19]. 


X shift = V-2in (r)cos (2nr 


a shift = tan(xr eer 


where r = random number between zero and one. 


ET 40 
Me 


These redistribution functions were tested by generating a million random 


numbers and plotting the resulting A shift (Figure 19 and 20). 
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Figure 19: Gaussian curve produced by redistribution function. 



















Lorentzian redistribution function 
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Figure 20: Lorentzian curve produced by redistribution function. 


Since collisional redistribution is dependent on whether an atom undergoes a 
collisional interaction, not all line-absorbed photons are collisionally redistributed. 
As a first order approximation, it can be assumed that the ratio of collisional de- 
excitation compared to total de-excitation is the same as the ratio of photons which 
undergoing collisional redistribution compared to just Doppler redistribution. 

It is with these two redistribution functions (eq 15 and 16) that re-emitted 
photons are redistributed in wavelength in SAM. The photons, once absorbed, are 
re-emitted with a new wavelength assigned randomly by a combination of these 
two functions. These functions are a major improvement over the simplified func- 


tions being used in SAM1_13. 
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IV. TESTING AND RESULTS 


Three methods of testing the validity of SAM were undertaken. The first 
method involved comparisons of theoretical mean free path lengths to computed 
mean free path lengths in a homogeneous atmosphere. The second method is com- 
parison of observed absorption line profiles to SAM-generated profiles. Finally, 
observational limb darkening is compared to SAM-generated limb darkening. To 
demonstrate the flexibility of SAM, stellar turbulence and stellar wind effects are 
applied to the solar Ti A6258.1 A line. 


A. MEAN FREE PATH 


In this test both the line absorption subroutine and the continuum absorption 
subroutines are tested. Although stellar atmospheres are not homogeneous, the 
homogeneous assumption allows mean free path calculations to be made. There- 
fore this is an important test for SAM. Mean free path calculations arc shown 
below. 


Treating the atmosphere as homogeneous and including only absorption 
processes of the continuum, the expected normal exponential decay equation is 


obtained. 


7 = exp(-x,pL) (17) 
0 


where: I = intensity at distance L through the medium 
Io = initial intensity, 
k, = continuum opacity per gram, and 
p = density. 

The average length (Lay. ) is; 











where x, and p are given by the photospheric model being used. 
For the mean free path of the line absorption process, the equation is; 


I = - 
i exp(-x,L). 
Therefore, 
1 
Lone = —— 
‘ave Ky, 
where 





f = oscillator strength, 

8, = statistical weight upper level, 

U = partition function of the atom, 

N; = population level of i state, 

6 = 5039.99/kinetic temperature, 

X = excitation energy of lower level (eV), 
Avp = Doppler shift, and 

H(u a) = Hjerting function. 


(19) 


(20) 


SAM generated L,,- values were produced by averaging 3000 path lengths 


through the homogeneous medium. Table 1 and 2 contains various input parame- 


ters as well as both calculated L,,, and SAM-generated L,,, . 


Mean free path 


comparisons for the continuum subroutine are done for four different continuum 


optical depths, while the line absorption subroutine are done for four different line 


optical depths. 
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TABLE 1 :CONTINUUM 
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TABLE 2 :LINE 
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A comparison of the results reveals the good agreement between the typical 
mean free paths calculated by SAM and those predicted by theory. The minor 
difference between the two values may by accounted for by the averaging of the 


random processes. 
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B. COMPARISON TO OBSERVATION 


The second method of validation is the comparison of the computer generated 
lines to observed lines. Line observations are taken from the National Solar Obser- 
vatory Atlas [Ref. 20). 


1. Convergence 


Figures 21, 22 and 23 show a time series from SAM. Since SAM is a 
Monte Carlo simulation, the more photons SAM generates, the higher the resulting 
signal to noise ratio in the SAM generated curve becomes. The trade off can be 
seen in these three figures. All SAM line profiles were generated on a Sun 
SPRAC1 workstation and CPU times stated are from this computer. Figure 21 
required only 2 minutes of CPU time and produced 10,000 photons. Figure 22 
required 20 minutes of CPU time and produced 100,000 photons. Figure 23 
required about 200 minutes of CPU time producing 1,000,000 photons. 
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Figure 21: 10,000 photons produced by SAM in 2 minutes 
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Figure 22: 100,000 photons produced by SAM in 20 minutes. 
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Figure 23: 1,000,000 photons produced by SAM in approximately 200 
minutes. 
Standard fluctuation calculation reveals that the fluctuation per bin is in 
agreement with random processes. 
(20) 


F = Nove 
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where: F = Fluctuation per bin, and 
Nave = Average number of photons per bin. 


2. Line synthesis 

One of the major applications of SAM is the fitting of stellar and atomic 
parameters. For this thesis, three solar line are synthesized. The portions of the 
solar spectrum being synthesized are shown in Figures 24, 25 and 26. The figures 
show two sets of lines. The dark line is the intensity of the solar flux on a scale 
between 0 and 100%. The lighter line is the more detailed scale running from 90 
to 100%. It is the darker line that SAM is modeling. The wavelengths given in 
these figures are in nanometers. These lines were chosen principally because of 
the availability of the two atomic constants; spontaneous emission coefficient (Ein- 
stein A), and quadratic stark interaction constant (C4). For comparison to SAM 
generated lines, the absorption lines near the absorption line being modeled have 
been replaced by a smoothed continuum. 

Observed portions of solar spectrum 
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Figure 24: Observed solar spectrum near Sr (4607.3 absorption line. 
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Figure 25: Observed solar spectrum near K (4044.1 absorption line. 
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Figure 26: Observed solar spectrum near Ti 46258.1 absorption line. 


Initial comparisons to estimated abundance, and C4 values are shown in 
Figures 27, 28, and 29. Notice that most of the fits lack the required line strength 
and line broadening. Abundance values are given as the log of the ratio of atom 


abundance to hydrogen abundance, with log H abundance = 12.0. 
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Comparison of observed lines to SAM-generated lines 
using initial C4 and abundance values 


4.0E+04 


3.0E+04 


flux 


2.0£+04 


1.0E+04 


—K sam generoted 
0-0 observations 


0.0 
4043.5 






4044.0 4044.5 
wavelength 


Figure 27: Observed K 14044.1 line to SAM-generated line with abun- 


dance = 5.6 and log C4 = -13.59. 
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Figure 28: Observed Sr 1 4607.3 line to SAM-generated line with abun- 


0.0 
4607.0 


dance = 2.7 and log C4 = -15.81. 
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Figure 29: Observed Ti A 6258.1 line to SAM-generated line with abun- 
dance = 4.78 and log C4 = -10.41. 


The initial values of abundance and C4 are adjusted until a good fit is 
obtained. Adjustment is appropriate as these values typically have a large uncer- 
tainties. Figures 30, 31, and 32 show SAM-generated fits with adjusted atomic 
values. Quality of the fits can be compared to that in Figure 3. 
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Comparison of observed lines to SAM-generated lines 


using adjusted C4 and abundance values 
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Figure 30: Observed K 4044.1 line to SAM-generated line with abun- 
dance = 5.81 and log C4 = -11.19. 
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Figure 31: Observed Sr (4607.3 line to SAM-generated line with abun- 


dance = 3.00 and log C4 = -9.81. 
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Figure 32: Observed Ti A6258.1 line to SAM-generated line with abun- 
dance = 4.80 and log C4 = -10.41. 


C. LIMB DARKENING 


Solar radiation intensity (ergs/cm”) decreases from the center of the solar 
disk to the limb. This limb darkening effect is caused by successive viewing of 
higher and higher photospheric layers, where temperature is lower and the source 
function is less intense. It is the limb darkening effect that allows temperature 
distributions and approximate source functions to be determined as a function of 
height in the stellar atmosphere. These are important parameters for stellar 
modeling. The understanding and successful simulation of this effect is therefore 
applicable and relevant. 

The SAM simulation retains two pieces of information on exiting photons: 
photon wavelength and exit angle from the surface normal. If SAM’s output is con- 
sidered spread uniformly over the solar disk then the angular binning of the exit- 
ing photons gives the exit location of the photons. But limb darkening measure- 
ments are taken from a thin spatial region from the center of the solar disk to its 


edge. These measurements in effect reduce the sun’s intensity into a 2 
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dimensional cross section, where the intensity is a function of angular position. 
But observations of limb darkening are made in terms of intensity (photons /cm 7”), 
therefore each angular region’s exiting photon flux must be divided by the projected 
surface area normal to the observer to obtain the necessary intensity. Figure 33 
shows a comparison of SAM-generated limb darkening for 5000A to observed 
curve. Agreement is good until the SAM-generated limb intensity becomes greater 


than observations around 70 degrees. 


SAM limb darkening compared to observations 
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Figure 33: SAM generated limb darkening at 5000A compared to 5000A ob- 


served curve. Deviation is still under investigation. 


Investigations into this effect seem to indicate it is related to the continuum 


opacity x, but earlier mean free path calculations indicate x, is calculated correctly. 
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D. STELLAR TURBULENCE 


Turbulence is the motion of the atmospheric gases on a scale that is large 
compared to atomic dimensions but small compared to the size of the star. Eddies, 
convection cells, and prominence-like structures are all forms of turbulence. These 
motions cause additional Doppler shifting to the photons and can greatly affect the 
line shape. 

The velocity Doppler shifting due to stellar turbulence occurs in SAM during 
the calculation of the line-absorption coefficient. SAM does not currently Doppler 
shift the photons due to stellar turbulence during re-emission. This would be fairly 
easy to correct if the nature of the turbulence is known. For the purpose of this 
thesis the re-emission correction is ignored. The broadening caused by increasing 
the stellar turbulence velocity to the Sr 14607.3 A is shown in Figure 35. 

SAM generated stellar turbulence 


flux 


m—a turbulence vel=4 Okm/sec 
iran SUrbUIONCS vale e- Cums aes 


4607.0 4607.5 4608.0 





wovelength 
Figure 35: SAM generated stellar turbulence for three different velocities. 
The top line is for a turbulence velocity of 6.0 km/sec, the middle line is for 
4.0 km/sec, and the bottom line is for 2.0 km/sec. Offsets are induced for 


viewing clarity. 
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E. STELLAR WINDS 


Stellar winds differ from turbulence in that the gases driven by convection or 
radiation obtain the necessary velocity to escape from the star. The amount of 
escaping matter can be quite large, typically 1 earth mass per year. The effects of 
this escaping matter on stellar lines is not well known. Although some general 
aspects of stellar wind effects on line profiles can be anticipated. The rapidly mov- 
ing matter will most certainly Doppler shift emitted photons. The Doppler shifting 
will only be in blue direction as the material which is moving away from the 
observer will be blocked by the star. Consequently, line profiles should contain 
increased flux in the shorter wavelengths. Many additional effects and the driving 
forces of these winds are not well known. This is one of the motivations behind the 
creation of SAM. 

In an effort to test the stellar wind routine in SAM, a 10 km/sec maximum 
stellar wind was imposed on the sun. Although this is clearly unrealistic, the goal 
was to reproduce some of the effects of stellar winds on a known spectrum. The 
wind speed varied linearly with altitude, reaching a maximum of 10 km/sec at the 
top of the photosphere. the equation for this variation is 


wind velocity = —2— maximum wind velocity , (22) 
zmax 


where z = distance out of the photosphere from reference depth (500km), 
zmax = distance to top of photosphere from reference depth. 


Figure 36 shows the effect of this wind on the earlier synthesized Ti absorption 
line. The large emission spike at the shorter wavelengths is the same effect 
observed in stars with large stellar winds [Ref. 21]. This is an exciting initial 
result. To emphasize the ease in which SAM was modified, only four additional 
lines of code were required to simulate movement of the outer photosphere. Two 


additional wind velocities were also run as shown in Figures 37 and 38. 
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Stellar wind effects on absorption lines 
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ma stellar wind mox vel = 10km/sec 
6258.0 6259.0 


wovelength 





Figure 36: SAM generated absorption line with 10 km/sec maximum stellar 





wind. 
Po 
2 
0.0! m—a stellor wind vel = 3 km/sec Tare white | 
6257.0 6258.0 6259.0 
wavelength 
Figure 37: SAM generated absorption line with 3 km/sec maximum stellar 
wind. “y 
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Figure 38: SAM generated absorption line with 20 km/sec maximum stellar 


“wind. 
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V. CONCLUSION 


The primary goal of this thesis has been the modification and testing of a new 
photon-following approach to stellar radiative transfer. The approach has several 
major advantages over standard stellar radiative transfer methods and has only 
recently been possible with the aid of fast computers. It was demonstrated that 
this technique accurately and efficiently reproduced observed solar absorption lines, 
including the effects of stellar rotation, stellar turbulence, collisional excitation, and 
collisional de-excitation. Stellar winds and limb darkening were reproduced 
although problems still remain in these areas. Nevertheless, the successful simula- 
tion of these effects is an improvement over previous methods. One of the major 
advantages of this approach to stellar radiative transfer is many vastly-different 


stellar phenomena can be modeled and incorporated into a single simulation. 


To demonstrate the usefulness of this technique, new abundance and C4 
values for three solar lines were determined. This technique can also be used to 
investigate almost all phenomena related to stellar atmospheres and radiative 
transfer (i.e., testing various model atmospheres, stellar wind theories, stellar tur- 


bulence, etc.) 


The results described in this thesis used fairly well-behaved stellar absorption 
lines. An extension of this thesis would be to generate true NLTE atomic popula- 
tion levels in the outer photosphere. If this were accomplished, the complex NLTE 
lines could be simulated. Some additional extension area are: accurate represen- 
tation of the Hjerting function when the Voigt a parameter exceeds unity, new 
redistribution functions that separate natural and collisional broadening, and 
improved velocity profiles for both stellar turbulence and stellar winds. Although 
these problems are non-trivial, there are no foreseeable road blocks for this new 


technique to simulate a large variety of complex stellar lines. 
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APPENDIX A. SOLAR ATMOSPHERIC PARAMETERS 


In charpter III section C spline fits to atmospheric parameters [Ref. 22] were 
replaced with interpolation tables. Figures 39 through 42 show a comparison 
between spline fit curves and interpolation values for temperature, density, opacity 
and gas pressure as each vary with depth. Incremental distance of 5km was 
chosen to allow the difference between the two curves to be seen. Better agree- 
ment is actually found in SAM because of the use of Ikm increments. Note an 
exception to the 5km increments is Figure 41 where an incremental distance of 
2km is chosen because of the reduced distances being plotted. 


11000 .K 


TEMP 


4000. K 
0. km DEPTH 600.00 km 


Figure 39: Temperature curve generated by spline fit compare to interpolation 
tables 
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Figure 40: Density (-4,) curve generated by spline fit compare to interpola- 
cm 


tion tables 
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Figure 41: Opacity a curve generated by spline fit compare to interpola- 


tion tables 
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Figure 42: Gas pressure re sone curve generated by spline fit compare to in- 
terpolation tables 
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